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A0535+26 is shown to display quantised IR excess flux states, which are inter- 
preted as the first observational verification of the resonant truncation scheme pro- 
posed by Okazaki and Negueruela (2001) for BeXRBs. The simultaneity of X-ray 
activity with transitions between these states strongly suggests a broad mechanism 
for outbursts, in which material lost from the disc during the reduction of truncation 
radius is accreted by the NS. Furthermore changes between states are shown to be 
governed by a ~1500 day period, probably due to precession of the Be disc, which 
profoundly dictates the global behaviour of the system. Such a framework appears to 
be applicable to BeXRBs in general. 



1 INTRODUCTION 



Roughly two thirds of High Mass X-ray Binaries possess a Be 
star as the primary, and are known as Be X-ray binaries or 
BeXRBs. A neutron star (NS) orbits the 09-B2 primary in 
10 — > 300days with pronounced eccentricity (e=0.1— >0.9). 
Whilst some systems display only weak persistent X-ray 
emission (X Per type), most additionally exhibit large tran- 
sient X-ray outbursts, sometimes close to their Eddington 
luminosity. Outbursts fall into two types (Stella et al. 1986): 

(I) Type I outbursts (10 36 - 10 37 erg s" 1 ) last a matter 
of days and occur close to or shortly after periastron. 

(II) Type II or 'giant' outbursts of ~ 10 38 erg s _1 occur 
between orbital phase and 0.5, lasting several tens of days. 

Nearly 3 decades after its discovery during a type II 
event, A0535+26 remains of considerable interest to those 
studying BeXRBs; due to its proximity high quality data 
can be obtained with ease, whilst in many ways it appears 
to be typical of its class. The 09.7IIIe (Steele et al. 1998) 
primary HDE 245770 displays typical variability in all as- 
pects of spectral, photometric and X-ray properties (Clark 
et al. 1998; Clark et al. 1999; Larionov et al. 2001). In addi- 
tion to type I outbursts which recur at 110-111 day intervals 
(presumed to be the orbital period P or b), 'giant' or type II 
outbursts have occurred less predictably in 1975, 1980, 1989 
and 1994. This paper is an attempt to draw together 15 years 
of multiwaveband variations to compose a coherent global 
picture of the behaviour of A0535+26. 

Integral to the analysis presented here is recent largely 
theoretical work regarding disc truncation in BeXRBs 
(Okazaki and Negueruela 2001; Okazaki et al. 2002). Central 
to the new picture is the abrupt truncation of the Be disc 
at a radius smaller than the periastron separation, where 



the orbital period is resonant with that of the neutron star 
and tidal interactions exceed viscous forces, the specific res- 
onance being dependant upon the system parameters and 
the viscosity a. The traditional picture of the NS passing 
through the Be disc is thus erroneous, explaining quiescent 
X-ray states in systems with robust discs. 

The mechanisms leading to accretion and X-ray out- 
bursts remain unknown except in systems of such high ec- 
centricity that inefficient truncation permits accretion at pe- 
riastron. In the case of 4U 0115+63 spectroscopy, photom- 
etry and X-ray observations have been combined to yield a 
picture in which quasi-cycles of 3-5 years govern the system 
(Negueruela and Okazaki 2001; Negueruela et al. 2001). In 
their interpretation, such cycles commence with disc-loss fol- 
lowed by reformation and expansion until the disc becomes 
truncated. The disc becomes unstable to warping/tilting, 
and precession causes alternating single peaked and shell 
line profiles. Type II outbursts then occur, presumably re- 
lated to the unstable disc configuration. 

This study was limited by the low S/N and sparsity 
of photometric and spectroscopic data for 4U 0115+63, at- 
tributable largely to the profound reddening. The analysis 
of A0535+26 presented here has the advantage of a 15 year 
dataset of spectra and IR photometry including many quasi- 
simultaneous spectral and photometric measurements and 
the opportunity to build upon the work that has been done 
in understanding 4U 0115+63. 

In Section 3 we present the first firm evidence for res- 
onant truncation of the Be disc with the discovery that the 
IR excess is quantised with relative fluxes matching theoret- 
ical predictions. In Section 4 we discuss a 1500 day cyclical 
'looping' behaviour seen in plots of WxboVS. tuk, as trun- 
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cation radius and emission line flux vary coherently but 90° 
out of phase. 

Section 5 describes how X-ray outbursts, in particular 
'giant' type II outbursts, occur at times of decreasing trunca- 
tion radius and at a specific cyclical phase. An association of 
'expelled' material with NS accretion is strongly suggested 
and a mechanism for X-ray activity suggested. A further 
point of interest uncovered in the dataset is the possible IR 
detection of the NS accretion disc previously inferred from 
X-ray QPOs (Finger et al. 1996) during the 1994 type II 
X-ray event. In Section 6 data obtained following the 1998 
disc-loss event are used in combination with the known disc 
size to determine the inclination, rotational parameter ui, 
and the circumstellar disc's radial velocity. 

Section 7 discusses the various periodicities observed 
in A0535+26, their probable origins and mutual associa- 
tions, in particular highly coherent photometric variations 
at the oft-detected 103 day period. Thus the phenomenology 
that has been described regarding 4U 0115+63 is applied to 
A0535+26 and some modifications suggested. The integral 
role of the warped, precessing disc invoked in 4U 0115+63 
is supported in the case of A0535+26, though in a less tran- 
sient role. We are thus led to a closer understanding of the 
global behaviour of the system. 
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2 OBSERVATIONS 
2.1 Photometry 

The Southampton/ Valencia collaboration has gathered in- 
frared observations at JHK from 1987 to 2001 with the 1.5-m 
Carlos Sanchez telescope (TCS), located at the Teide Ob- 
servatory in Tenerife, Spain. The data were reduced follow- 
ing the procedure described by Manfroid (1993). Instrumen- 
tal values were transformed to the TCS standard system 
(Alonso, Arribas & Martinez 1998). The data acquired for 
A0535+26 are plotted in Figure 1. 

For the histograms in Figure 2, multiple observations 
spanning several days from individual observing runs are 
represented by single data points, which are thus separated 
by at least 10 days to prevent statistical bias arising from 
unrepresentative data sampling. Typical photometric errors 
arc ±0.05. 

Photometric vtljhk measurements at epoch 17/12/1997 
have been extracted from the 2MASS catalogue at IRSA. 

The my photometric data presented by Lyuty and 
Zaitseva (2000), and available online as a VizieR catalogue, 
are used here and also plotted in Figure 1. Errors are ±0.005. 



2.2 Spectroscopy 

Spectra of A0535+26 including Ha and He I A6678 have been 
obtained since 1987 at numerous observatories. Observations 
not described in Clark et al. (1998) or Haigh et al. (1999) 
are listed in Table 1. The JKT, INT and WHT data have 
been obtained through the ING service programme or the 
INC Archive. 

Reduction took place using STARLINK packages, with 
spectral analysis performed in DIPSO. 

Twelve additional Ha equivalent width (W\H a ) obser- 
vations between Jan 1997 and Mar 1998 are taken from 



7.5 : 

8.0 : 



E 0-2 - 

o.o : 

-15 ! 

E : 
& -a! 

; 
x -1 - 



2447000 2448000 2449000 2450000 2451000 2452000 

Figure 1. Observational parameters of A0535+26 from 1987- 
2001. The upper panel shows detected X-ray activity and the 
cycles described in Section 4. 



Piccioni (1998). Two further values from 11-12/11/1998 are 
taken from Giovannelli et al. (1999), with one each from 
6 Nov 1998 (Lyuty and Zaitseva 2000) and 29 Jan 2000 
(Negueruela, private comm.). 
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Table 1. Unpublished spectroscopic observations of HDE 245770. 



Date 


Telescope 


Instrument 


Details 




13/12/1991 


INT 


IDS 


EEV5, 500 camera, 


R1200R 


8-10/3/1993 


SAAO 1.9m 


Spectrograph 


RPCS 




12/2/1995 


JKT 


AGB 


EEV5, AGBX 




4/4/1996 


SAAO 1.9m 


Spectrograph 


RPCS 




28/10/1997 


JKT 


AGB 


TEK4, AGBX 




6/2/1998 


SAAO 1.9m 


Spectrograph 


SITe, 1200 1mm- 1 




24/12/1998 


WHT 


ISIS blue arm 


EEV12, R600B 




27/3/1999 


WHT 


ISIS red arm 

lUllJ 1 l^-VA Cil Ill 


TEK2, R600R 




21/9/1999 


INT 


IDS 


TEK5, 235 camera, 


R1200R 


20/10/1999 


WHT 


UES 


SITel 




22/3/2000 


INT 


IDS 


TEK5, 235 camera, 


R1200R 


10/4/2000 


INT 


IDS 


TEK5, 500 camera, 


R1200Y 


15/8/2000 


INT 


IDS 


TEK5, 500 camera, 


R1200Y 


16/10/2000 


Mt. Skinakas 1.3m 


Spectrograph 


1302r 1mm" 1 +80fim + ISA 608 


5/2/2001 


INT 


IDS 


TEK5, 500 camera, 


R1200Y 


11/11/2001 


SAAO 1.9m 


Spectrograph 


SITe, 1200 1mm- 1 
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Figure 2. Histograms of A0535+26 decretion disc magnitudes 
m disc ; n j ^ k and colour index (J-K). Bin widths are 0.2 for 
md jHK anc ^ 0-1 for m (j S ^ K y Flux states A and B are labelled on 
the mj^ sc histogram. 



3 IR EXCESS QUANTISATION AND 
TRUNCATION 

In order to study only the circumstellar decretion disc flux 
the constant stellar component m* is subtracted - using stel- 
lar magnitudes m} = 8.6, m* H = 8.5 and m* K — 8.5 - and the 
resultant magnitude subsequently indicated with the super- 



script disc. Note rnj l ^ c K are relatively insensitive to errors 
in m* JHK . 

Figure 2 shows histograms of the vrij^ K ij-k) photome- 
try obtained from 1987 to 2001. These infrared fluxes display 
a clear bimodality as the system alternates between faint 
and bright flux states denoted A and B respectively and sep- 
arated by Am* 1 " magnitudes. The contrasting constancy of 
m {j-K) snows the dominance of optically thick isothermal 
emission (Dougherty et al. 1994), suggesting that changes 
in disc area are responsible, rather than optical depth or 
temperature. 

Such a quantized luminosity distribution is a natural 
consequence of the resonant truncation hypothesis proposed 
by Okazaki and Negueruela (2001). Here tidal torques, which 
operate most strongly in orbits resonant with that of the 
NS, truncate the disc at radii meeting the condition Pns ■ 
Ptmnc = n : 1 where n is an integer. They find that for 
viscosities of a = 0.013 — 0.38 the disc is truncated at the 
6:1 resonance, a = 0.038-0.11 at 5:1 and for a = 0.11-0.40 
the disc is supported out to 4:1. 

Assuming isothermal optically thick emission, the quan- 
tity Am* sc is a direct measure of the relative areas of the 
disc in these two states. Knowing the orbital period of the 
disc's outer edge at a given resonance, and assuming Ke- 
plerian rotation and a non edge-on inclination enables the 
relative disc radii and thus Am dlsc to be predicted inde- 
pendently of stellar mass M». Thus a change in truncation 
resonance produces a Am dlsc plotted in Figure 3, dependant 
on the size of any central cavity, or separation of the inner 
disc radius from the photosphere. 

Rivinius et al. (2001) found that in the typical isolated 
Be star r\ Cen, mass injection events or 'outbursts' lowered 
the disc inner radius to 1.2i?„ , but at other times a consistent 
3.1 — > 3.3i?* was observed. There is no reason to believe that 
the same mechanisms should not lead to similar cavities in 
A0535+26. Curves for different values of r inner assume i?* = 
14.7.R© and M» = 26M (Vacca et al. 1996) based upon an 
09.7IIIe classification; note the star's volume produces a 
'cavity' of ri nner — 1R* . Note that resonances are with the 
orbital beat frequency of Pb ea t = 103 days rather than the 
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Figure 3. Modelling of disc flux ratios Am* 81 for changes be- 
tween disc states truncated at ni:l and n2:l resonances. Emission 
is from an optically thick isothermal disc with central cavity of 
firmer = 1 , 4, 4.3, 5R* , where R* = 1.02 X 10 10 m (Vacca et al. 
1996). Assuming the identities of states with resonances discussed 
in the text, the A symbols show the observed values of Am* sc . 



orbital period of P or b — HI days for reasons discussed in 
Section 7, though this has a minimal effect. 

Of the three IR photometric wavebands K has the 
largest optical depth and most closely approximates the op- 
tically thick assumption. Photospheric emission is also min- 
imized in this band, so the Am* sc data are used hereon, 
though all wavebands yield similar results. Table 2 and Fig- 
ure 2 yield a value for Am*' c of 0.60; consulting Figure 3 
this suggests a central cavity Tinner / R* — 4.3 taking states 
A and B to be truncated at 6:1 and 5:1 respectively. Whilst 
the various assumptions make this only an approximate re- 
sult, M» and 7?, are sufficiently well determined to yield 
maximum errors in r inner I R* of ~ 30% requiring a cen- 
tral cavity to be present. Lack of isothermality will cause 
Tinner to be underestimated if as expected Tdisc decreases 
radially. Departures from optical thickness will primarily af- 
fect the outer disc and will have the same effect. This re- 
sult is in encouraging agreement with the central cavities 
observed in r\ Cen, whereas using radii for dwarves yields 
Tinner / R* — 7.5. This result therefore supports the giant 
(III) classification. 

IR photometry (Persi et al. 1979), and the my data 
presented by Lyuty and Zaitseva (2000), suggest that the 
disc was in an even brighter state (C) from the early 1970s 
until 1980 (see Table 2). Optically thick emission in the IR 
bands requires this to have been a physically larger disc. 
As Okazaki and Negueruela (2001) note, truncation at 4:1 
places the outer edge of the disc marginally outside the 
Be star Roche lobe at periastron, and thus represents the 
largest stable disc state. The identity of state C with this 
4:1 truncation, and thus states B and A with 5:1 and 6:1, 
is appealing as it explains the prodigious X-ray activity of 



Table 2. Observational parameters and ranges for the two ob- 
served flux states from 1987 to 2001, the additional state C pro- 
posed from the literature and disc-loss. l:Using ms of Lyuty 
et al. (2000) and m [B _ v) = 0.55, 2: 11-13/12/1976, Persi et 
al.(1979), 3: Haigh et al. (1999), 4: Voikhanskaya (1980), Baratta 
ct al.(1978). 



State C B A Discloss 



w XHa 


-17 ±3 4 


-12 ±4 


-10 ± 3 


~ 2.5 


w XHeI 


? 


-0.2 ±0.3 


-0.8 ± 0.4 


0.3 ±0.2 


m v 


~ 8.85 iO.l 1 


9.10 ±0.1 


9.30 ±0.1 


9.45 3 


mj 


~ 7.34 2 


7.60 ±0.1 


8.00 ± 0.15 


8.49 3 


m H 


~ 7.02 2 


7.35 ±0.1 


7.75 ± 0.15 


8.39 3 


m K 


~ 6.71 2 


7.10 ±0.1 


7.50 ±0.15 


8.34 3 


mdisc 

m K 


~ 6.94 


7.45 ±0.1 


8.05 ± 0.15 


> 10 



the system during the 1970s, periastric Roche lobe overflow 
alone enabling accretion and X-ray activity (Okazaki and 
Negueruela 2001). Identifying state C with a 3:1 truncation 
places much of the disc outside the Be star Roche lobe at 
periastron preventing stability, and requires an unrealistic 
viscosity (Okazaki et al. 2002). 

Figure 3 shows that a change in truncation between 5:1 
and 4:1 (state B^C) will produce a Am* sc « 0.56 and thus 
predicts rnj^ sc ~ 6.89 for state C, in excellent agreement 
with the isolated measurement of mf sc = 6.94 ±0.03 in the 
literature (Table 2). 

It should be noted here that the strongly shepherded 
structure of a truncated disc leads to highly stable contin- 
uum fluxes - this is discussed in Section 7. 



4 CYCLICAL BEHAVIOUR 

Figure 4a plots ttik vs. W\Ha between 1987 and 2001. For 
each measurement of W\Ha the closest photometric vhk 
measurement - always within 30 days and mostly within 10 
- is taken. Flux states A and B are clearly visible as parallel 
bands of points. 

A0535±26 can be seen to follow a clockwise loop with a 
mean period of 4.36 ±0.5 years (1590 ± 180 days) except for 
those points around the disc-loss episode (Haigh et al. 1999) 
at lower left. Following Negueruela and Okazaki (2001) 
who observed similar cyclical behaviour in the BeXRB 
4U 0115±63, these looping cycles in A0535±26 are defined 
as commencing following disc-loss, or at the state A— >B tran- 
sition which occurs at the same phase. In terms of the flux 
states described above, each cycle consists of approximately 
equal periods first in state B followed by state A. During 
state B, W\Ha gradually increases before decreasing follow- 
ing the B — > A transition. Figure 5 shows photometry and 
X-ray activity during cycle 2, which has the best observa- 
tional coverage and appears to be typical. 

In contrast to W\Ha which varies through a similar 
range in both states, Figure 4b shows W\Hei oscillating be- 
tween two loci in phase with the WA_f/ Q variations. With the 
state B locus at W\H e i — —0.15 and W\uei ~ —0.60 in state 
B, the emission line flux ratio of states A:B~ 4, demonstrat- 
ing that the smaller changes in continuum strength (~ 1.3, 
determined from tub) cannot be responsible. Allowing for in- 
filling of the underlying photospheric absorption He I A6678 
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Figure 4. Path of m K vs. W X Ha (a) and VK AHe/ A6678 (b). Flux 
states A and B are labelled. In a), movement around the loop is 
clockwise. 

has a ~ 25% greater flux in the fainter state A. This is inter- 
preted as arising from enhanced densities at the small radii 
where this line is primarily produced (Stee et al. 1998), due 
to the more strongly confined disc in this state. At the same 
time, the wings of He I A6678 arising in the inner disc dis- 
play emission to systematically higher velocities. Thus these 
diagrams enable the flux state to be determined either from 
spectra or photometry. 

Due to the patchy nature of the dataset, a holistic anal- 
ysis of optical and IR photometric and spectral data was re- 
quired to reconstruct the timing and phenomena associated 
with the ongoing and previous 3 cycles, which are marked at 
the top of Figure 1. Whilst the ~ 1400 day photometric com- 
ponent of these cycles has been previously noted in period 
searches (Clark et al. 1999; Larionov et al. 2001), the co- 
herent variation of other observational parameters has not. 
These observations demonstrate that the long period vari- 
ations observed in A0535+26 are not arbitrary but rather 
linked to a 1400-1500 day cyclical variation, perhaps pre- 
cessional, in the system, the nature of which is discussed in 
Section 7. 



5 TEMPORAL ASSOCIATION WITH X-RAY 
OUTBURSTS 

When the disc outer radius makes a quantum leap inwards 
(state C— >B or B— >A) the material present between those 



radii must be relocated elsewhere in the system. It is sug- 
gestive to note that such a state change was observed in mx 
between JD 2449344 (state B) and 2449401 (state A), while 
a type II X-ray outburst was detected by BATSE (Finger 
et al. 1996) from 2449380 to 2449432. This simultaneity is 
well seen in Figure 5, and strongly suggests that the giant 
outburst was the accretion of material previously resident 
between the 5:1 and 6:1 resonances. The greater duration 
of variability in ray relative to tuk is interpreted as being 
due to the greater optical depth in rriK desensitising the 
photometry to small variations in disc structure as trunca- 
tion progressed. The entire 1993-4 series of X-ray outbursts, 
types I and II, all occurred during the somewhat lengthier 
truncation process as seen in mv, and presumably corre- 
spond to material accreted following removal from the outer 
disc. Importantly, truncation would appear to be the cause 
of the X-radiation and not vice-versa, as truncation com- 
mences before X-ray emission. 

In 1989, precisely one cycle earlier during the state B — ► 
A transition of cycle 1, a type II followed by a type I X-ray 
outburst occurred in the system. Whilst observations are 
sparser, this series of events appears extremely similar to 
the 1993-4 events described above, and we therefore suggest 
that the same processes were responsible. 

Unlike cycles 1 and 2, the type II outburst in Oct. 1980 
occurred during state C with my = 8.9 ±0.04 from Jan 1980 
until Dec 1980. No state change occurred until some months 
after the event. However this can be understood, as the 4:1 
resonance hypothesised to bound state C is marginally out- 
side the Be star Roche lobe (see Section 3) enabling mass 
transfer without a change in truncation radius (Okazaki and 
Negueruela 2001). This may indeed be the mechanism be- 
hind many of the X-ray outbursts during the bright photo- 
metric state of the 1970s, including the 1975 discovery type 
II event. As has been noted (Lyuty et al. 1989) even the qui- 
escent X-ray flux of A0535+26 was elevated from 1975-1980; 
the identification of this epoch with a 4:1 truncation (state 
C) permits mass transfer through the inner Lagrangian point 
at periastron even in the absence of truncation changes. The 
absence of accurate IR photometry from this period prevents 
further analysis at the present. 

Clark et al. (1999) in discussing the long term photo- 
metric variations of A0535+26, noted that X-ray outbursts 
seemed to occur after a period in which the optical light 
curve was fading. Clark suggested that this represented an 
episode in which disc material was lost radially, with X-ray 
emisison triggered by the interaction of this material with 
the NS. However, this suggestion was rejected in favour of 
the traditional 'shell ejection from the primary' scheme. 

The question remains of why some transitions lead to 
X-ray activity while others do not. A possible factor is that 
though disc truncation may occur at a specific phase in the 
1500 day precession cycle, the precession angle changes by 
26° between each periastron passage so system geometries 
do not recur precisely. To take a simplistic viewpoint, it may 
be that whether the NS lies precisely in the disc plane at the 
time of critical precession phase determines whether a type 
II outburst occurs or not. The 1996 truncation occurred very 
early in the cycle, and no X-ray activity was seen. 

As with the model of Negueruela et al. (2001) for 
4U 0115+62, this model of X-ray outbursts is effectively in- 
dependant of Be star activity, and eliminates the unwieldy 
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Figure 5. Photometry throughout cycle 2 including the 1993-4 
X-ray outbursts. The duration of the BATSE 20-100 keV detec- 
tion during the type II event is shown as a short line. 



'ejections' of material from the Be star frequently invoked 
as directly fuelling X-ray outbursts (de Loore et al. 1984). 
Such connections between long-term (months-years) photo- 
metric changes and X-ray events has not received a great 
deal of attention in the literature, with much work expect- 
ing to observe short optical flares associated with enhanced 
stellar mass- loss prior to an outburst. 

It should be noted here that whilst photometry is un- 
available, spectroscopy from 2 Jan 2003 showing W\H a =- 
13.5A and W XHe i A6678=-0.3A (Laycock, private comm.) 
confirms that the system is in the latter stages of state 
B (see Table 2). The imminent transition to state A, ex- 
pected mid 2003, could produce the first X-ray outburst in 
the A0535+26 system since 1994. 

5.1 IR excess colour changes 

The presence of an accretion disc around the neutron star 
in A0535+26 has been debated in the literature (Giovan- 
nelli and Ziolkowski 1990; Lyuty and Zaitseva 2000). While 
no evidence has been found in the optical/IR to date, X- 
ray observations have provided unequivocal evidence during 
the 1994 giant outburst. Finger et al. (1996) detected 20- 
lOOmHz QPOs and large P from JD 2449380 to 2449432 
which they concluded are due to an accretion disc around 
the NS. Motch (1991) and Li (1997) found that an accretion 
disc is only present during such giant outbursts, therefore 
the data obtained at this time presents an excellent oppor- 
tunity to look for emission from any such disc. 

Figure 5 plots my, niK and the IR colour index m( J _ /f ) 
for the duration of cycle 2. In particular it shows an anoma- 
lous transient decline in m(, / _ /f ), from JD 2449401 to 



2449440, simultaneous with the 1994 state B— >A change, 
associated type II X-ray outburst and known NS accretion 
disc. Gaps in the data span 57(182) days before(after). Such 
a run of points at a low rri(j-K) relative to the typical value 
of 0.52 ± 0.05 is unique in the 14 year dataset. 

This large change in near IR fluxes cannot be due to the 
underlying B star, and the required increase in T disc while 
the overall disc flux is decreasing makes little sense in terms 
of normal models of disc emission. Another possibility is that 
being simultaneous with a state change, transient emission 
from residual low density material at radii between the new 
and old truncation radii may alter the global IR properties. 
However, such emission, assuming T < Tdisc, is expected 
to have redder colours than the optically thick disc bulk 
(Dougherty et al. 1994) and so cannot be responsible. 

An additional transient source of emission with a rela- 
tively blue colour index is thus suggested, the additional flux 
compensated for by the simultaneous reduction in Be disc 
emission associated with the state B — > A change. If so an 
intriguing possibility is that this is emission from the known 
accretion disc around the NS. Limits can be found for the 
additional flux in m,/, the primary uncertainty being the Be 
disc flux. Based purely upon the colours, even if none of the 
additional flux is present in vtlk the 0.1 magnitude colour 
shortfall requires 10% additional flux in m,j. Alternatively 
if we assume the Be star to have faded fully to a typical 
early state A luminosity then the total system magnitude 
excluding any accretion disc is m^ e = 7.45 ±0.10 (Figure 4) 
giving rrij e = 7.97 ±0.11 (in excellent agreement with the 
value in Table 2) assuming a typical rnfj_ K - j of 0.52 ± 0.05. 
The observed mj = 7.80 ± 0.07 therefore yields a maximum 
excess of 0.17±0.13 in mj. Thus a 10-20% excess in mj with 
a bluer m,( j_ K ) than the dominant circumstellar emission is 
found. 

Approximating the radius of the accretion disc to the 
NS Roche lobe radius at periastron and adopting the canon- 
ical Mns = 1.4M© and M, = 26M (Vacca et al. 1996) 
yields r L(NS) = 3.45 x 10 lo m (Eggleton 1983). Such a disc 
has 22% of the emitting area of the Be disc in state A, which 
in turn emits 44% of the J band flux. If we further assume 
optically thick emission and temperature Tad = Tsedisc an 
additional 9.7% flux in mj is predicted. Therefore the es- 
timate of a 10-20% excess in m.j, combined with the bluer 
m (.j-K), strongly suggests a significantly higher Tad and, 
possibly, a smaller emitting region. 

The simultaneity of this unique set of photometric data 
within the 14 year dataset with an independent signature 
of an accretion disc in the system is noteworthy. The tran- 
sient existence of an accretion disc in A0535±26 is thus sup- 
ported. 



6 DISC EXPANSION AND INCLINATION 

Spectra obtained following the 1998 disc-loss episode pro- 
vide a valuable glimpse of a reforming disc. Figure 6 plots 
W\Ha, the separation of the peaks of the Ha profile AV peafe , 
and the derived disc radius r oute r- The fit to the data in the 
lower panel shows that r oute r can be modelled as increas- 
ing linearly at 824ms _1 with a break around JD 2451200. 
During this break, disc expansion halts and W\H a remains 
constant; this presumably represents a hiatus in disc refor- 
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Figure 6. Ha line parameters following the 1998 disc-loss 
episode. Lower panel assumes i = 36.1° and M* = 26Mq. 



mation. Such a highly subsonic expansion velocity is in good 
agreement with standard models of viscous Keplerian discs. 

Following disc-loss in 4U 0115+63, disc expansion con- 
tinues until truncation becomes effective (Negueruela and 
Okazaki 2001). In A0535+26 this can be seen to have oc- 
curred near JD 2451650; ma shows state A (Figure 1) and 
thus truncation at the 6:1 resonance (though note a bright- 
ening to state B and 5:1 shortly afterwards). 

Assuming Ha emission to the outer edge of a rotat- 
ing disc, the Doppler displacement v obs of each peak of the 
resultant two peaked line profile represents the line of sight 
component of the Keplerian velocity at r ou ter, whilst the res- 
onant truncation paradigm informs us of the orbital period 
p at this radius. Stee and Bittar (2001) found that for typi- 
cal isolated Be stars the Ha emitting region was bounded at 
18-R*. As confirmed below, for all reasonable values of M* 
the disc is fully within this region in both states A and B. 

The product of the orbital velocity and period yields 

T outer 

2Txr outer = ^^rP (1) 
aim 



whilst for a particle in a circular orbit 



sin 2 i , . 
r = — = — GM, 



(2) 



Eliminating r then gives a measure of i as a function of 



Vobs 



2-kGM, 



1/3 



(3) 



which is plotted in Figure 7 for all reasonable stellar 
masses. The case of the state A truncation at JD 2451650 
described above is used, and Ha AV pea k=2v bs is mea- 
sured at 288.9 ± 2kms~ 1 , in excellent agreement with the 
275 ± 20 km s^ 1 measured during the 1994 state A. For the 
reasons discussed in Section 7 the apparent NS orbital pe- 
riod relative to the disc is Pteat = 103 days, though this 
almost no effect on the result. Section 3 established that 
state A is probably truncated at the 6:1 resonance giving 
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Figure 7. i vs. M*/Mq for two identifications of the truncation 
resonance in state A. Note 6:1 is favoured by theory and obser- 
vation, discussed in Section 3. 



p = 17.17 days, though i vs. M, is also plotted assuming a 
5:1 commensurability yielding inclinations ~ 3° larger. 

Uncertainties stem largely from the uncertainty in M* . 
The standard spectral classification 09.7IIIe requires a mass 
of 26±3M Q (Vaccaet al. 1996), yielding 36. 1°±1. 7; the less- 
favoured dwarf (V) classification (M* = 22 ± 3Mq) yields 
38.6° ± 2.5. These are at the low end of previous published 
results (Janot-Pacheco et al.1987). 

Returning to Equation 1, r outer = (5.79±0.2) x 10 10 m at 
a 6:1 truncation and (6.68±0.25) x 10 m at 5:1. Depending 
upon luminosity class and using R* from Vacca et al. (1996) 
5:1 truncation corresponds to r ou ter = 6.79i?«(11.9-R*) for 
luminosity class HI(V) verifying the assumption of Ha emis- 
sion over the disc's full extent. This result clearly under- 
scores the extent of truncation in BeXRBs relative to the 
discs of isolated Be stars. 

Of course the rotating disc is not necessarily constrained 
to lie precisely in the equatorial plane of the star; Section 
7 strongly suggests that it does not. However the constancy 
of the measured line velocities through the cycles requires 
that the tilt cannot be large (Ss 10°) and so whilst this 
assumption should be born in mind the inclination should 
be accurate to approximately this tolerance. 



6.1 Determining uj 

The disc-loss episode of 1998 provided a rare chance to 
obtain spectra free from circumstellar emission. Following 
Mazzali et al. (1996) blue end He I lines were used to measure 
the stellar v sini due to their intrinsically narrow FWHMa, 
using photospheric spectra obtained shortly after disc-loss at 
the WHT. The preferred line, narrow with minimal emission, 
at 4713.20A yields a vsini of 251 ± 5kms 1 , whilst the in- 
trinsically broader A4471,A4922 lines give 254.2 kms -1 and 
261.3 km s~ respectively. Removing instrumental broaden- 
ing, determined from arc lines to be 120 km s -1 , yields a 
vsini of 220.4kms~ 1 at A4713, and 224 (232) kms -1 for 
A4471 (A4922). 
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Thus a value of v sini= 225 ±10 km s" 1 is adopted. This 
compares with the published value of 250 ± 15kms _1 by 
Wang and Gies (1998), who did not have the advantage of 
spectra uncontaminated by emission. A giant(dwarf) classi- 
fication and thus i = 38.6° ± 2.5(i = 36.1° ± 1.7) implies an 
equatorial rotational velocity of 361 ± 25(382 ± 20) kms -1 . 

The quantity uj gives the ratio of equatorial rotational 
velocity to 'critical' velocity, that of a particle in a circu- 
lar orbit at the stellar equator. Porter (1996) found that 
Be shell stars show ui — 0.67 ± 0.04 and provided over- 
whelming evidence that this is true of Be stars in general. 
Depending upon a giant/dwarf classification, the masses 
and radii of Vacca et al. (1996) yield values for v cr u of 
582kms~ 1 (719kms~ 1 ) for luminosity class III(V), giving 
w = = 0.62 ± 0.05(0.53 ± 0.03). This is consistent with 
a giant but excludes dwarves, supporting the 09.7IIIe clas- 
sification. 



7 PERIODICITIES AND A GLOBAL MODEL 

Whilst many period searches have been undertaken using 
photometry of A0535+26 (see Larionov et al.(2001), hence- 
forth L2001, for a review), only 103 and 1400 day periods 
repeatedly emerge at high significance whilst the orbital pe- 
riod is most notable by its absence. The 1400 day period is a 
manifestation of the cyclical photometric states described in 
this work, but the origin of the 103 day period is less clear. 

Whilst the strongly shepherded structure of a truncated 
disc leads to highly stable IR photometric fluxes in all states, 
in the most strongly confined state (A) even mv demon- 
strates remarkably little random variation. This photomet- 
ric stability is highly conducive to period searching. Thus 
54 ray data spanning 468 days during state A of cycle 2 
(1994-5) were searched for periods in the range 90-120 days 
using the PDM software (Stellingwerf 1978). The resultant 
periodogram in Figure 8 again finds a strong 103.1±0.1 day 
periodicity and nothing at P or b- Folding using this period, 
the 102.83 ± 0.15 day period of L2001 and two published 
orbital periods, reveals the remarkable lightcurves in Figure 
9. 

While perturbations at the 103 day period result in disc 
flux changes by a factor of ~ 2 in my, random fluctuations 
are < 0.03 magnitudes, and may in fact be due to obser- 
vational error. Thus the luminosity, and thus presumably 
the disc structure, is determined entirely by the NS orbital 
phase. This may also be the cause of the larger range of tuk 
in state A than state B (Figure 4. 

It is clearly of considerable relevance that the beat pe- 
riod of the orbital with the 1400 day period is 103 days. 
As L2001 suggest this is strongly suggestive of a 1400 day 
precessional period in one of the major flux producing sys- 
tem components. Adopting a range of mv = 9.26 —> 9.36 
from Figure 9, and using a photospheric mv = 9.45 ± 0.02 
from the 1998 disc-loss episode (Figure 1 and Haigh et al. 
(1999)) gives a range of ~ 2 for the Pbeat fluctuations in disc 
emission. This variation is in severe conflict with the limit 
of a few percent to the total mv flux found by Lyuty and 
Zaitseva (2000) for the possible contribution of an accretion 
disc. Furthermore the apparently transient nature of the NS 
accretion disc (Section 5.1) argues against its involvement. 

Therefore the precession must occur in the Be decretion 
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Figure 8. PDM periodogram of 54 my data spanning 468 days 
from state A of cycle 2 (1994-5). Resolution is 0.05 days, using 
6 PDM bins. Note strong detection at 103.1 ± 0.1 days and non- 
detection at the ~ 111 day orbital period. Data from Lyuty and 
Zaitseva (2000). 
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Figure 9. mv photometry from state A of Cycle 2 (1994-5) 
folded at the detected beat period and two published orbital 
ephemcrides. Data from Lyuty and Zaitseva (2000). 



disc, in a retrograde sense to the NS orbital motion in order 
to yield a NS-disc perturbation period shorter than P or b- 
Indeed for a disc rotating in the same sense as the NS such 
retrograde precession is required (Wijers and Pringlel999). 
Simulations by Okazaki et al. (2002) provide an example of 
how the disc's emission can then be modulated in mv at 
the 103 day period. However because their calculations of 
phase resolved disc flux assumed optically thin ff/bf emis- 
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sion, the probable optical thickness in mv during state A 
will alter the variations quantitatively but retain the peri- 
odicity. In addition to confirming lack of variability at P or b, 
the apparent perturbations occurring at 103 days suggests 
that NS-disc interactions occur at this shorter period rather 
than Port- For this reason a 103 day period is used as the 
'apparent' orbital period Pb ea t in the flux modelling of Fig- 
ure 3. 

Section 4 showed how A0535+26's observation parame- 
ters, most notably the truncation radius, have in recent years 
varied in well defined cycles of mean duration 1590 ± 180 
days. Using the detected 103.1±0.1 day period as a beat pe- 
riod of P orb = 111.0 ±0.5 with a third period gives 1450 ±90 
days, consistent with the mean cycle duration; this equality 
is strongly suggestive of a common origin. The cycle pe- 
riod derived from measurements of Pbeat and P OT b is prob- 
ably more accurate than simply timing photometric state 
changes; the larger value found from cycle durations appears 
largely due to delayed increases in truncation radius as the 
Be disc declined around the 1996 minimum. 

The 1590 day cycles of truncation radius are explicable 
as arising from geometric variations in the strength of tidal 
torques as a function of precession angle. A possible clue 
to such continuous tidal variation is the trend from square- 
wave to quasi-sinusoidal variation at ~1500 days towards 
shorter continuum wavelengths (compare mv to uik in Fig- 
ure 1 through cycles 1 and 2). Whilst tuk represents emis- 
sion from the surface of a well defined disc my, in general 
only partially optically thick, is sensitive to disc density and 
therefore reflects the continuously changing effect of tidal 
forces upon mass distribution in the disc. 

A constant precession of the decretion disc - as opposed 
to transient precession prior to X-ray activity - is therefore 
the main suggested alteration to the model of 4U 0115+63 
put forward by Negueruela et al. (2001). 

The identification of specific truncation resonances al- 
lows constraints to be placed on the viscosity parameter a 
by comparison with the simulations and theoretical work of 
Okazaki and Negueruela (2001). The geometrically varying 
tidal forces will most closely match these theoretical predic- 
tions when the Be disc and NS orbit are most closely copla- 
nar. At this time such forces will also be strongest and thus 
maximum truncation will occur; therefore state A, truncated 
probably at 6:1, occurs in the most coplanar phase. Thus for 
stability to occur at this truncation radius suggests the con- 
straint a < 0.013 - 0.038 (Okazaki and Negueruela 2001). 

7.1 A general model for BeXRBs? 

The processes described in this work are fundamentally 
an extension of the theories put forward regarding the 
BeXRB 4U 0115+63 (Negueruela et al. 2001; Negueruela 
and Okazaki 2001). However, photometry of the BeXRB 
4U 1145-619 presented by Stevens et al. (1997) clearly shows 
X-ray outbursts occurring during a fading in mv after a 
spell at a constant luminosity, an anticorrelation apparently 
identical to that described for A0535+26 in Section 5. 

Approximately 4 yr periods of disc growth/dissipation 
strikingly similar to those reported here for A0535+262 have 
been reported for several other BeXRBs. Reig et al. (2001) 
present photometry and spectroscopy of LS992/RX J0812.4- 
3114 showing a ~ 4.3 year cycle; following disc-loss IR ex- 



cess rapidly builds until the onset of X-ray activity, with 
IR excess fading slightly as Balmer emission peaks - this is 
identical to the activity cycles of A0535+26 seen in Figures 
1 and 4. 

The dataset on which the model for 4U 0115+63 (Sec- 
tion 1, Negueruela et al. (2001)) is based also shows near- 
identical variation of Balmer lines, IR excess and X-ray ac- 
tivity over 3-5 year cycles. 'Quasi-cyclic Be star envelope 
variations' have been invoked to explain the 4 year modu- 
lation in the amplitude of the periodic radio outbursts in 
LSI+61°303 (Gregory et al. 1989), while Hummel (1998) re- 
ported shell-Be-shell transitions in 7 Cas at 4 ± 0.4 years 
which he attributed to a precessing circumstellar disc. Clark 
et al. (2001) also show that X Per has exhibited disc- loss and 
disc 'low states' separated by intervals of 5 and 6 years. 

If disc precession is the root cause of all these cycles, 
the similarity of cyclical period in systems with disparate 
orbital periods and inclinations is intriguing. While the disc 
is probably ousted from the equatorial plane by NS interac- 
tions, attribution of the precession to the asymetric potential 
of the oblate OB star is an appealing explanation, as the pri- 
maries of BeXRBs are known to occupy only a small range in 
u> (Porter 1996) and M, (Negueruela 1998) and thus should 
possess almost identical gravitational properties. 

Behaviour similar to that observed in A0535+26, LS992 
and 4U0115+63 is probably restricted to systems of mod- 
erate eccentricity causing the disc to be truncated at a low 
resonance (i.e. 4:1) so that changes between adjacent reso- 
nances cause observable flux changes. The degree of copla- 
narity of NS orbit and Be equator will clearly be an influen- 
tial parameter profoundly affecting the nature of the tidal 
interactions. Orientation on the sky is also doubtless impor- 
tant for understanding different systems; it seems likely that 
the Be disc of 4U 0115+63 intercepts the line-of-sight at a 
specific cyclical phase, causing transient shell lines in that 
object and not in A0535+26. Clearly further multiwaveband 
studies of other BeXRB systems will demonstrate how gen- 
eral or indeed accurate this model is, particularly as differ- 
ences emerge. 



8 SUMMARY 

The primary findings of this study are the detection of quan- 
tised disc fluxes strongly supporting the resonant truncation 
paradigm, association of X-ray outbursts with reductions in 
truncation resonance, and the role of a ~1500 day Be disc 
precession period in governing this system behaviour. Op- 
tical observations have revealed variability at the 'effective' 
Port, relative to the precessing Be disc, and represent Be 
disc perturbations by the NS. The standard luminosity clas- 
sification, III, is supported. 

This greater understanding of the A0535+26 system ap- 
pears to be applicable to BeXRBs in general, perhaps allow- 
ing deeper insights into Be discs and quasi-Keplerian discs 
in general. 
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